% Plot the power at M2 (12.421 hours) of the observatory data on to the
% predicted power map


S = dir('/home/mnair/projects/tides/*_2010.mat');
load '/home/mnair/projects/tides/M2_Abs_Colormap' map1;
open('/home/mnair/projects/tides/M2_Bz_Abs_0_720_1440.fig');
hold on


for nstn = 1:length (S);
    
  eval(['load /home/mnair/projects/tides/' S(nstn).name ]);

  L = 1./(F*3600) >= 12.419 & 1./(F*3600) <= 12.4240;  
  
  if any(L)
  m2_amp = mean(sqrt(P(L,1)/(en-st)));
  
  % Find color to plot
  
  cidx = floor(63*(m2_amp)/(10)+1);
    
    if cidx > 64,
        cidx = 64;
    elseif cidx < 1,
        cidx = 1;
    end;

    plot(longitude, 90-latitude, '.', 'color',map1(cidx,:),'MarkerSize',50);
  end;
    
end;